Neutron-proton interaction in rare-earth nuclei: Role of tensor force 
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We investigate the role of the tensor force in the description of doubly odd deformed nuclei within 
the framework of the particle-rotor model. We study the rare- earth nuclei 174 Lu, 180 Ta, 182 Ta, and 
188 Re using a finite-range interaction, with and without tensor terms. Attention is focused on the 
lowest K = and K = 1 bands, where the effects of the residual neutron-proton interaction are 
particularly evident. Comparison of the calculated results with experimental data evidences the 
importance of the tensor-force effects. 
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1 I. INTRODUCTION 

x> _ 

It has long been known that spectroscopy of doubly odd deformed nuclei evidences phenomena which can be directly 
associated with the interaction between the unpaired neutron and proton. The two most important effects of this 
kind are known as Gallagher-Moszkowski (GM) splitting [1] and Newby (N) shift [2]. 

Recently, use of large multi-detector 7-ray arrays with high analyzing power has largely extended the amount of 
experimental data on this class of nuclei. For instance, currently available data in the rare-earth region allow the 
empirical determination of 137 GM splittings and 36 N shifts for 25 doubly odd nuclei [3,4]. These are to be compared 
with the 50 GM splittings and 19 N shifts known in the mid 1970s, as reported in the extensive review article by 
Boisson et al. [5]. The new experimental data have also evidenced the occurrence of other phenomena which may 

■ provide further information on the neutron-proton (n-p) interaction. Of particular interest is the odd-even staggering 
in K ^ bands, the most important mechanism responsible for it being the direct Coriolis coupling with one or more 
N-shifted K = bands [6] . Also to be mentioned is the signature inversion phenomenon, which in the last few years 
has attracted much attention leading to several experimental and theoretical studies. Various attempts have been 
made to understand this phenomenon in term of different mechanisms, but its interpretation still remains an open 
question [7]. 

In a previous paper [8] (hereafter referred to as I) we have studied the doubly odd deformed nucleus 1 Lu by 
performing a complete Coriolis band-mixing calculation within the framework of the particle-rotor model. We focused 
^ I attention on the lowest K v = 0~ and K v — 1 + bands, where the effects of the residual n-p interaction are particularly 

■ evident producing a significant N shift in the former and an odd-even staggering in the latter. The aim of our study 
was in fact to investigate the role of the tensor force, which has long been a controversial matter. A detailed discussion 
on this point can be found in I. In that paper we have used for the n-p interaction both a finite-range force with 
a Gaussian radial shape and a zero-range interaction. To completely explore the role of the tensor force, we have 
performed two different calculations with the Gaussian potential, with and without the tensor terms, respectively. 
The main result of these calculations was that the tensor force is absolutely essential for the description of the above 
mentioned bands in 176 Lu. In I we have also confirmed the relevance of exchange forces, which was already emphasized 
in Refs. [9,10]. 

Here, we extend our calculations to other doubly odd deformed nuclei in the rare-earth region, focusing attention 
on N-shifted and odd-even staggered bands in 174 Lu, 180 Ta, 182 Ta, and 188 Re (some preliminary results have already 
been reported in Refs. [11-13]). The main motivation for this work is to verify if the conclusions based on the study of 
176 Lu are confirmed when considering a larger set of experimental data. In fact, the role of the tensor force in doubly 
odd deformed nuclei seems not yet completely recognized, as evidenced by the large number of recent calculations 
making use of zero-range forces (see for example Ref. [14] and references therein). 

The outline of the paper is as follows. In Sec. II we give a brief description of the model and some details of our 
calculations. Our results are presented and compared with the experimental data in Sec. III. In Sec. IV we draw the 
conclusions of our study. 
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II. OUTLINE OF THE MODEL AND CALCULATIONS 



As mentioned in the Introduction, our calculations are performed within the framework of the particle-rotor model, 
in which the unpaired neutron and proton are strongly coupled to an axially symmetric core and interact through a 
residual effective interaction. This is a well known model and a detailed description can be found, for instance, in 
Refs. [5,15]. For the sake of completeness, we give in the following the most relevant formulas. The total Hamiltonian 
is written as 

H = H + Hjipc + H ppc + V np . (1) 

The term H includes the rotational energy of the whole system, the deformed, axially symmetric field for the neutron 
and proton, and the one-body intrinsic contribution from the rotational degrees of freedom. It reads 

H, = ^(I 2 - /f) + H n + H p + - f n3 ) + (f p f p3 )]. (2) 

The two terms -Hrpc and H ppc in Eq. (1) stand for the Coriolis coupling and the coupling of particle degrees of 
freedom through the rotational motion, respectively. Their explicit expressions are 

Hkvc = -^{I + J-+I-J + ), (3) 



and 



h 2 



H pP c — 2j-tinjp + 3n.jp)- (4) 

In Eqs. (2)-(4), J is the moment of inertia of the core while I and J = j p + j„ are the total and intrinsic angular 
momentum operators, respectively. I3 and J3 are their projections on the intrinsic symmetry axis and, owing to the 
axial symmetry, are represented by the same quantum number K . 
The effective n-p interaction in (1) has the general form 

V np — V(r)[uo + u\ + u 2 Pm + u 3 Pm + VtS 12 + VtmPmS 12 }, (5) 

where the notation is just the same as that adopted in Ref. [5]. 

As basis states we use the eigenfunctions of Hq, which are written as a symmetrized product of rotational functions 
D I MK and intrinsic wave functions \v n Sl n ) \v p Sl p ) 



f 21 + l\ ^ 

\v n n n V p n p I M K) = (-j^2~) [ D MK\ u n£ln)\Vp£lp) 



+ {-) I+K D I M _ K \ Vn n n )\ Vp n p )\. (6) 

Here SI is the the quantum number corresponding to j'3 and v stands for all the additional quantum numbers necessary 
to completely specify the states. The state \vQ) is the time-reversal partner of \vQ). The quantum number K has 
two possible values 

K± = \\{\ n \±\n p \\, (7) 

corresponding to parallel or anti-parallel coupling of Q n and Q p . 

Each intrinsic state gives then rise to a rotational band, but the Hamiltonian (1) produces an admixture of different 
bands. The explicit expressions of the matrix elements of the total Hamiltonian can be found in Ref. [15], where it 
can be seen that while the Coriolis interaction has only AK = ±1 matrix elements, the two terms V np and -ff ppc give 
rise to diagonal and non-diagonal contributions. In particular, the latter has diagonal matrix elements different from 
zero only for K = bands with \Q n \ — \Sl p \ — 1/2. We would like to stress that all diagonal and non-diagonal terms 
of Hamiltonian (1) arc explicitly taken into account in our calculations. 

As regards the single-particle Hamiltonians H n and H p , they have been generated by a standard Nilsson potential 
as defined, for instance, in [16]. In our calculations for the four nuclei 174 Lu, 180 Ta, 182 Ta, and 188 Re the parameters 
fi and k in this potential have been fixed by using the mass-dependent formulas of Ref. [17], while the harmonic 
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oscillator parameter v has been chosen according to the expression v = muj/Ti fm~ 2 . The deformation parameter (3% 
has been deduced for each doubly odd nucleus from the even-even neighbor. In Table I, we report the values of /?2 
together with those of the rotational parameter, h 2 /2J , used in our calculations. For each doubly odd nucleus this 
latter quantity has been deduced from low-lying bands with a pure rotational character. The adopted single-particle 
schemes for the neutron and proton are listed in Table II for the four considered nuclei, each scheme being essentially 
derived from the experimental spectra of the two neighboring odd-mass nuclei. 

As regards the residual n-p interaction (5), we have used a finite-range force with a radial dependence V(r) of the 
Gaussian form 

V(r)=exp(-r 2 /r 2 ). (8) 

As mentioned in the Introduction, the main aim of this study is to assess the role of the tensor force. To this 
end, we have performed two different calculations with the Gaussian potential, with and without the tensor terms, 
respectively. In both cases, for all the parameters of the force we have adopted the values determined by Boisson et 
al. [5], the only exception being u 0} which does not contribute either to the GM spittings or to the N shifts. These 
values, as well as that of uq, are reported in I, where the reasons of our choice are also discussed. 

Before closing this Section it is worth mentioning that we have also performed calculations making use of a 5 
interaction with a strength of the spin-spin parameter varied over a large range of values. We have found that no 
value of this strength gives a satisfactory description of the N shifts and odd-even staggerings in any of the considered 
nuclei, thus confirming our results for 176 Lu (see I). As a consequence, here we will not compare with experiment the 
results obtained by using a S force, but make only some comments to point out the inadequacy of this force to explain 
the above mentioned effects. 



III. RESULTS AND COMPARISON WITH EXPERIMENT 

We report here the results of our calculations for 174 Lu, 180 Ta, 182 Ta, and 188 Re, which are all well deformed nuclei. 
The choice of these nuclei is motivated by the fact that enough experimental information relevant to our study is 
available for them. In fact, in the experimental spectra of 174 Lu, 182 Ta, and 188 Re K — bands, which arc not 
strongly perturbed but have a sizeable N shift, have been recognized, while in 180 Ta and 188 Re there are K = 1 bands 
which exhibit a significant odd-even staggering. We have focused attention on those data which are directly related to 
the n-p interaction, so as to give a sound answer to the question regarding the role of the tensor term. In fact, aside 
from the residual n-p interaction, also the terms H ppc and i?R,pc of the Hamiltonian (1) may give a contribution to 
the odd-even shift (see Rcf. [6]) and a complex interplay of all these interactions may occur in bands characterized 
by a strong admixture of different intrinsic wave functions. In these cases, even taking explicitly into account all 
contributions coming from Hamiltonian (1), it would be a very difficult task to disentangle the genuine effects of the 
n-p interaction. 

In Figs. 1-3 the experimental spectra of the K = bands [4] are reported and compared with the calculated ones, 
which are obtained making use of the two different interactions mentioned in Sec. II. The K + = + p|[404]n|[633] 
band of 174 Lu starting at 281 keV excitation energy is shown in Fig. 1. This band, recognized in several other odd-odd 
nuclei, has been observed in 174 Lu up to I n = 9 + . We sec that a very good agreement with experiment is obtained 
for the Gaussian force with tensor terms. This is not the case for the calculations performed with the pure central 
interaction. From Fig. 1 it appears that the two calculations yield almost the same level spacings for states with even 
and odd angular momenta separately. The main difference between the two spectra resides in the energy displacement 
of levels of odd I relative to those of even /, i.e. the N shift. In first order perturbation theory the N shift is just 
defined as the matrix clement (i> p Q, p , v n U, n \V np \v P ^l P ^ v n Qn)- In Table III we report the values of this matrix clement 
obtained with the central plus tensor force together with the individual contributions. For the above mentioned K = 
band in 174 Lu we see that the magnitude of the matrix element is almost entirely determined by the tensor force, 
the central part giving only 1 keV contribution. This result is in line with the early findings of Rcfs. [5,18], where 
it was evidenced the importance of the contribution of tensor force to the N shift for bands based on antiparallel or 
parallel intrinsic spins. For the latter the effects of the tensor part are predominant, as expected from the study of 
Ref. [2], where an analysis in the asymptotic limit was performed. In this context, we should also mention that when 
using a 5 force, as defined for instance in I, there is no reasonable way to reproduce the odd-even shift of this band. 
In fact, taking a negative value for the strength v\ of the spin-spin term, the odd-even shift has the opposite sign 
as compared to the experimental one, its magnitude increasing almost linearly with v\. As an example, when using 
v\ =-0.9 MeV, which is a reasonable value to reproduce the GM splitting in the rare earth region [3,19], the matrix 
element (v p Q. p , v n Q, n \V np \v p Q, p , v n Qn) becomes 22 keV. The fact that a 6 force produces a shift with the wrong sign 



3 



while a pure central force gives a very small contribution, but with the right sign, makes evident the importance of 
the exchange forces. 

In Figs. 2 and 3 we compare the experimental and calculated spectra for two lowest K = bands in 182 Ta and 
188 Re, respectively. The first is the K + = CP p|[404]n|[503] band with a bandhead energy at 558 keV while the 
second is the K + = + p|[514]n|[505] band starting at 208 keV. Unfortunately, only low-spin members of these two 
bands have been identified, but, as for 174 Lu, an excellent agreement is obtained in both cases when a Gaussian force 
with tensor terms is used. From Figs. 2 and 3 we see that the right level order is obtained while the discrepancies 
in the energies are less than 20 keV for all the states. As regards the calculations with a pure central force, the 
situation is quite similar to that discussed for 174 Lu. In fact, the disagreement between the experimental data and the 
results of these calculations resides in the relative displacement in the energy of the states with odd and even angular 
momentum, which, as mentioned before, is essentially due to n-p interaction. In Table III we report the values of the 
relevant matrix elements. We see again that the pure central force gives only a small contribution, the magnitude of 
these matrix elements being essentially determined by the tensor force. 

To conclude this discussion, we have shown that the tensor force is quite able to explain the odd-even shift in 
K = bands of 174 Lu, 182 Ta and 188 Re. We now compare in Table IV the empirical GM splittings [3] for these three 
nuclei with those derived from the calculated bandhead energies. In this Table we only report the results obtained 
using the Gaussian plus tensor force. This is to show how this force, parametrized as in Ref. [5], leads to a very good 
description of the three nuclei 174 Lu, 182 Ta and 188 Re on the whole. In fact, we see from Table IV that not only is the 
sign correctly reproduced in all cases, also the quantitative agreement is quite satisfactory. The discrepancy between 
the calculated and experimental values is always less than 50 keV, except in one case where it turns out to be 101 
keV. In this case, however, the empirical GM splitting is affected by an error of 57 keV. 

Let us now come to the odd-even staggering effect. In Fig. 4 the experimental [20,21] ratio [E(I) — E(I — 1)]/2J 
for the K n = l + p|[404]n|[624] ground-state band in 180 Ta is plotted vs I and compared with the calculated ones. 
We see that this band exhibits a rather large odd-even staggering, the even I states being energetically favored with 
respect to the odd ones. This behavior is satisfactorily reproduced by the calculation including the tensor force while 
with the pure central Gaussian force the staggering is almost nonexistent. It should be pointed out that the staggering 
in this band may be traced to direct Coriolis coupling with the K 7 * — + p|[404]n|[633] band. In fact, in both our 
calculations we have found that the wave functions of the states of the if" = 1 + band contain significant components 
(up to 30%) of states of the above K w = + band. The fact that only the calculation including the tensor terms gives 
the right staggering implies therefore that only this force is able to produce a sizeable N shift for this K n = + band. 
Since this band has not been recognized in 180 Ta, a direct comparison is not possible at the present time. It should 
be noted, however, that the K * — + band in 174 Lu (see Fig. 1) has just to the same intrinsic n-p configuration. As 
we have already pointed out above, only with the Gaussian plus tensor force the energy of the even / states of this 
band is decreased so as to bring the calculated spectrum in good agreement with the experimental one. 

In Fig. 5 the experimental staggering observed in the = l~p|[402]n§[512] ground-state band of 188 Re [4] is 
compared with the calculated ones. Although this band does not exhibit a large odd-even staggering, we see that only 
with the Gaussian plus tensor force the calculated behavior reproduces the experimental one. Also in this case the 
wave functions of the states of the K n = 1~ band obtained from both calculations contain appreciable components (5 
- 10 %) of states with K n = 0~p§[402]n|[503]. The difference in the two calculations may therefore be traced back 

to the different values of the matrix element (p|[402]n|[503]|V r Tl p|p|[402]n|[503]). In fact, as can be seen from Table 
III, only the matrix element of the Gaussian plus tensor force turns out to be adequate to produce a sizeable N shift 
for the = 0~ band. Unfortunately, this band has not been observed in 188 Re or in other rare-earth nuclei. 

IV. CONCLUDING REMARKS 

In this paper, we have presented the results of a study of the well deformed doubly odd nuclei 174 Lu, 180 Ta, 182 Ta, 
and 188 Re within the framework of the particle-rotor model. We have performed complete Coriolis band-mixing 
calculations, focusing attention on the lowest K = and K = 1 bands, where the effects of the neutron-proton 
interaction are particularly evident. As regards this interaction, we have used in our calculations a Gaussian force 
with and without tensor terms adopting in both cases the values of the parameters originally determined by Boisson 
et al. [5]. 

This study represents an extension of our previous work on 176 Lu aimed at obtaining further insight into the role 
of the tensor force in the description of doubly odd deformed nuclei. The results presented here, as those reported 
in I, show that the tensor force is an essential ingredient, if not unique in some cases, to explain the odd-even shift 
in K = bands. Also the staggering in K ^ bands, whose basic mechanism is the Coriolis mixing with N-shifted 
K = bands, is very well accounted for by our calculations with tensor terms. 



4 



It should be mentioned that in more complex situations, where there is a strong band mixing, other mechanisms 
may play a significant role and only the right interplay between them can explain the observed effects. This is, for 
instance, the case of the so called doubly decoupled bands, namely of bands based on high- j unique-parity shell-model 
states. Also in this context, however, the tensor force should not be neglected, as is done in most of the existing 
studies to date. The use of an incomplete neutron-proton interaction may in fact lead to an incorrect evaluation of 
the weight of other mechanisms. 
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FIG. 1. Experimental and calculated spectra of the lowest K* = + band in 174 Lu. The theoretical spectra have been 
obtained by using (a) a central plus tensor force with a Gaussian radial shape, (b) a Gaussian central force. 



FIG. 2. Same as Fig. 1, but for the lowest K" = 0" band in 182 Ta. 

FIG. 3. Same as Fig. 1, but for the lowest K* = 0+ band in 188 Re. 

FIG. 4. Experimental and calculated odd-even staggering of the K n = 1 + ground-state band in 180 Ta. Solid circles correspond 
to experimental data. The theoretical results are represented by open circles (Gaussian central plus tensor force), and squares 
(Gaussian central force). 



FIG. 5. Same as Fig. 4, but for the K n = 1 ground-state band in 188 Rc. 
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TABLE I. Values of the parameters /3 2 and |j (keV). 



Nucleus 


h 


n 2 

2.7 


1V4 Lu 


0.331 


12 


180 Ta 


0.278 


14 


182 Ta 


0.274 


14 


188 Re 


0.226 


18 



TABLE II. Energies (in keV) of the proton and neutron single-particle states. 



Nucleus 



Configuration 
Proton 



E 



Configuration 
Neutron 



Nucleus 



Configuration 
Proton 



E 



Configuration 
Neutron 



2 

9- 
2 



IF 

| "[514] 
± [510] 
11512] 
i"[521] 
11512] 
| [503] 

I + [633] 
| + r, 



174 



Lu 



[404] 
[541] 
[402] 
[514] 
[411] 
[532] 





208 
369 
437 
449 
903 



5 
2 

2 + 

2 
7 
2 



11512] 
[633] 
[514] 
11521] 

11510] 
[642] 



5 + 
2 




389 
400 
431 
1053 
1207 



J Ta 



[404] 
[514] 
§1402] 
i + [411] 
11541] 





20 
250 
546 
700 



[624] 



-[642] 




200 
395 
518 
647 
785 
812 
1506 
2000 



2 Ta 



7 - 
2 


"[404] 





i - 

2 


[510] 





9 - 
2 


■[514] 


5 


3 - 
2 


[512] 


220 


5" 

2 


h [402] 


350 


11 - 
2 


H [615] 


300 


1 - 
2 


-[411] 


477 


94 

2 


[624] 


550 








7 - 
2 


[503] 


670 








7 - 
2 


[514] 


1100 



'Re 



[402] 



11514] 
i + [411] 
I + [404] 





176 

625 
748 



11512] 
i [510] 

7 
2 
9 
2 



_ [503] 
_ [505] 
§-[503] 





161 
310 
340 
640 



TABLE III. Values of the matrix elements (v p Q p , v n Q, n \Vnp\u p Q, p , v n Q, n ) (in keV) obtained using 
a central plus tensor force with a Gaussian radial shape. (C) and (T) represent the contributions 
of the central and tensor force, respectively. 





Configurat 


ion 








Nucleus 


Proton 


Neutron 


(C) 


(T> 


(Vn P ) 


m Lu 


| + [404] 


| + [633] 


-1 


-26 


-27 


182 Ta 


| + [404] 


11503] 


3 


19 


22 


188 Re 


11514] 


| [505] 


-4 


-60 


-64 


188 Re 


| + [402] 


§-[503] 





47 


47 
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TABLE IV. Experimental and calculated GM splittings (keV). 



Configuration 



Nucleus 


Proton 


Neutron 




K 71 


Expt. 


Calc. 


1 74 t 


3 [404] 


f [512] 
J [633] 
§[521] 
f [512] 
§[512] 


1 




-114.9 ± 1.3 


-124 




I [404] 


0+ 


7+ 


-64.2 ± 10.0 


-82 




5 [404] 


3~ 


AT 


79.5 ± 5.0 


89 




|[541] 
| [402] 


2 + 


3 + 


-149.8 


-123 




0" 


5~ 


130.2 ± 2.1 


155 




§[514] 


§[512] 


2 + 


7 + 


135.3 ± 3.1 


109 


i! "Ta 


I [404] 
I [404] 


±[510] 
§[512] 
| [503] 


3" 


AT 


-92.6 ± 8.3 


-118 




2" 


5" 


128.0 ± 0.7 


142 




|[404] 


0" 


7" 


-122.0 ±3.1 


-143 




§[514] 


±[510] 


4+ 


5+ 


147.6 ± 0.6 


105 




|[514] 


§[512] 


3+ 


6+ 


-111.1 ± 12.4 


-123 




§[514] 


#[615] 


1" 


10" 


147.7 ±0.6 


122 




§[402] 


3 [510] 
§[512] 


T 


3" 


301.2 


296 


188 Re 


|[402] 


1" 


4" 


-138.7 ± 1.4 


-117 




§[402] 


§[510] 


2" 


3" 


93.5 ±1.1 


100 




§[402] 


| [503] 


1" 


6" 


209.4 ±0.2 


205 




3 [4H] 


3 [512] 


1" 


2" 


197.2 ± 57.3 


96 



8 



800 


9+ 


782 


9+ 


735 






700 










9+ 


689 


600 


8+ 


562 
















8+ 


539 


8+ 


531 


500 


7+ 


501 










E(keV) 






7+ 


449 






400 










7+ 


401 


300 


6+ 
- 5+ 


313 
295 


6+ 


294 


6+ 


289 








5+ 


248 






200 


3+ 


164 






5+ 


199 




4+ 


1 A n 


4+ 


129 


4+ 


126 


100 


1 _ 


102 


3+ 


- 123 












1+ 


63 


3+ 


73 




2+ 


39 


2+ 


36 


2+ 


35 





0+ 





0+ 





i — 


11 





Expt. 


(a) 




(b) 















5" 


423 


400 


5" 
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